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A contribution (^^ to the curvature perturbation will be generated during the waterfall that ends 
hybrid inflation, that may be signiflcant on small scales. In particular, it may lead to excessive 
black hole formation. We here consider standard hybrid inflation, where the tachyonic mass of the 
waterfall field is much bigger than the Hubble parameter. We calculate (^^ in the simplest case, and 
see why earlier calculations of are incorrect. 



I. INTRODUCTION 

The primordial curvature perturbation C is one of the most important features of the 
early universe. On each scale it is relevant until the era of horizon entry at the epoch 
k = aH, and provides the principle (perhaps the only) initial condition for the subsequent 
evolution of all other perturbations. At the beginning of the known history of the universe, 
when the temperature is around 1 MeV, observation has established the existence of ( at 
wavenumbers between k ~ aoHo (corresponding to the size of the observable universe) and 
k ~ e^^aoHo. On these 'cosmological scales' its spectrum V(^{k) is almost scale-invariant 
with value "P^ (5 x 10-^)^. 

The primordial curvature perturbation presumably exists also on much shorter scales and 
at much earlier times. Although not directly observable, this can have a significant effect on 
the evolution of the early universe. The most dramatic possibility is black hole formation; 
assuming that C, is nearly gaussian, cosmological constraints on the abundance of primordial 
black holes requires Vq ^ 10~^ at horizon entry, on practically all scales [ij. Another 
important effect may occur if there is an era of matter domination in the early universe. 
A significant primordial curvature perturbation may cause the formation of gravitationally 
bound objects, with possibly dramatic effects 

The initial condition for the observable universe is presumably set by an early era of 
inflation, which at the classical level makes the universe absolutely homogeneous. The 
curvature perturbation must then originate as a quantum fluctuation. 

To explain the nearly flat spectrum, the observed curvature perturbation is generally as- 
sumed to originate from the vacuum fluctuation of one or more light scalar fields during 
almost exponential inflation. The curvature perturbation generated in this case can exist 
down to the scale k = e'^°aoHo, the horizon scale at the end of inflation. With a typical 
cosmology and inflation scale, A^o — 60. The spectrum of the curvature perturbation gener- 
ated by light field perturbations might increase strongly on small scales, allowing significant 
black hole formation [3]. 

Heavy scalar fields can acquire a classical perturbation through what is called preheating. 
Standard preheating (sj invokes an interaction g^cjP'x^ , where x is the heavy field and cj) is an 
oscillating field that is initially homogeneous. It begins with an era, during which remains 
homogeneous so that the field equation for Xk is 

Xk = -(m2(t) + A:2)xk, (1) 



* A preliminary version of this paper appeared as larXiv: 1005. 24611 
We use the standard cosmology notation. The comoving wavenumber fc of a Fourier component defines a 
comoving scale 1/fc so that smaller scales have bigger k, but k itself is loosely referred to as the scale. The 
physical wavenumber is k/a(t) and the Hubble parameter is H = a/a, while oq and Hq are evaluated at 
the present epoch. During inflation scales 'leave the horizon' when aH exceeds k and after inflation they 
they enter it when aH falls below k. 

A scalar field a with canonically kinetic term is said to be light during inflation if |rn^| H^ , where 
= d^V/dcr^ and V is the scalar fleld potential, and heavy if |m^| ^ H^. A vector field perturbation 
might also contribute to ^ 0|. 
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with rn^it) — rn? + g'^(fP'{t) and m? the mass-squared of x- If 'rn?(t) is dominated by the 
second term, the solution can grow exponentiahy, converting the vacuum fluctuation of 
Xk to a classical perturbation. The contribution to the curvature perturbation generated 
during this era is calculated in |6|| assuming the inflation potential 1^ cx 0^, and in 7] for 
the smooth/mutated hybrid inflation P potential. 

In this paper we consider instead tachyonic preheating, that ends hybrid inflation. Here 
the interaction is the same as before, but the waterfall field x has a negative (tachyonic) mass- 
squared. If (j) fa-Us below a critical value during inflation, m?{4>) becomes negative and Xk 
can grow exponentially, again converting its vacuum fluctuation to a classical perturbation. 
To facilitate the calculation we make some fairly restrictive assumptions about the waterfall, 
which will be relaxed in future papers. 

The layout of the paper is as follows. In Section |TT] we recall the definition and properties 
of C- In Section Hill we recall the basics of hybrid inflation. In Section HVl we lay down our 
assumptions and describe the evolution of waterfall field. In Section |V] we define the region 
of parameter space in which our assumptions are consistent. In Section IVTl we calculate the 
pressure and energy density of the waterfall field. In Section lVIII we justify our assumptions. 
In Section rVIIII we calculate C,^, the contribution of the waterfall field perturbation to the 
curvature perturbation, by integrating the perturbation theory expression for C,. In Section 
HXlwe review the 5N approach and show that it reproduces the previously-calculated result 
for C^. In Section |X] we see that previous calculations of C,^ are incorrect. In Section IXll we 
conclude, pointing to the need for a better understanding of the ultra-violet cutoff in the 
context of cosmological scalar field calculations. 



II. PRIMORDIAL CURVATURE PERTURBATION C 
A. Cosmological perturbations 

Let us recall some basic concepts, described for instance in [9|. Given a function /(x, 
in our Universe one can define the perturbation 5f: 

f{^,t)^f{t) + Sf{^,t), (2) 

where f(t) refers to some Robertson- Walker (unperturbed or background) universe whose 
line element is 

ds^ = -dt"^ + a^{t)5^jdx*dx' = a^irf) {-drf + S^dx^dx^) . (3) 

In the background universe the energy density p and pressure p are related by the continuity 
equation 

p = -3H{p + p), (4) 

where H = d/a. Also, Einstein gravity is assumed corresponding to the Friedmann equation 
p = 3M^ij2 with Mp = (87rG)-i/2 = 2 x 10^^ GeV. 

The coordinates (x, t) label points in the background universe, and also points in our 
Universe. The latter labeling (gauge) defines slices (fixed t) and threads (fixed x) in our 
Universe. Let us denote 5/(x, t) in some chosen gauge by f/(x, t). We will need the first-order 
gauge transformation for the case that / is specified by a single number. Going to a new 
coordinate system with time coordinate t(t,x), it is 

9-9 = -f{i~t) (5) 
= fSt, (6) 

where 6t is the time shift going from a slice of uniform t, to one of uniform t with the same 
numerical value. 

We will also need to describe the statistical properties of g at fixed t, which are simplest 
in Fourier space; 

gk(t)= J d3xe-*-^g(x,i). (7) 
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When considering the statistical properties one ignores the spatial average (zero mode), 
which can be absorbed into the background unless / corresponds to anisotropy. 

Considering an ensemble of universes, we assume that the the statistical properties are 
invariant under translations (statistical homogeneity) and rotations (statistical isotropy) 
By virtue of the former, ensemble averages can be replaced by spatial averages within the 
articular realization of the ensemble that corresponds to our Universe (ergodic theorem) 

Statistical homogeneity implies (gt) = 0. The two-point correlator defines the spectrum: 

(.gk.9p) = {2Ti)H\k + Y>){'2ii^/k'')V,(k). (8) 

An equivalent quantity, also called the spectrum, is Pg = / k^)V g . 

A gaussian perturbation is defined as one with no correlation between the Xk, except that 
implied by the reality condition x_k = Xk- only correlators are the two-point correlator, 
and the disconnected 2" point correlators starting with 

(5kigk2 5k3 5k4) = (ffkiSka) (ffkaSki) + permutations. (9) 
Non-gaussianity of g is specified by its bispectrum Bg, trispcctrum Tg etc., defined by 

(5ki5k..gk3) = (2^)353(ki+k2 + k3)i?g, (10) 
(5ki5k2<?k3 5k4>c = {2^f5^{ki+'k2 + 'k3 + '^i)Tg, (11) 

etc., where the subscript c denotes the connected part of the correlator that has to be added 
to the disconnected part. 

The mean-square perturbation {g^) can be taken as the spatial average of <7^(x), and is 
given by 

1 /• r°° rlk 

^ (^ y ^'^^^(^^ = L T^^(^) 

After smoothing g on a scale L this becomes^^ 

{g') = I ^Vg{k). (13) 







For a gaussian perturbation, the probability distribution of (7(x) is gaussian so that (g") 
vanishes for odd n and is equal to (n— l)!!((g^)"^^ for even n. For a non-gaussian perturbation 

= (2^)-6y"d3^id3fc2B<,(ki,k2), 

(/)^ ^ (/) - i{g^f = (2^)-9 j £k,dVhTgik,,k2, kg), (14) 

etc.. The skew, kurtosis etc. of the probability distribution are Sg = {g^) / {g^)^^^ , Kg = 

(5^)0/(5^)^ etc.. If they are all <C 1 g can be regarded as almost gaussian. 

We will also be interested in the spectrum of g^. Using Eq. © with the convolution 
theorem 

(3')k = 7:7^ / d'fc'.gk'.9k'-k, (15) 



(27r)3 



one finds [IC 



Pg^ik) = / d'k'Pgik')Pg{\k~k'\). (16) 



A function is said to be smooth on the scale L if its Fourier components are neghgible in the regime 
k > L^^ . Smoothing means that we remove Fourier components in this regime. 
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B. Primordial curvature perturbation ^ 



In the gauge with the comoving threading and the shcing of uniform energy density p, the 
spatial metric defines Q non-perturbatively through 

g,,(x,t) ^ a\^,t)(^^(-A , (17) 
\ / ij 

C = In {a{x,t)/a{t)) = Sin a. (18) 

In this expression Tr/i = 0, so that the last factor has unit determinant and a(x,t) is the 
local scale factor such that a comoving volume element is proportional to a^. 

For this definition to be useful, we need to smooth the metric and the stress-energy tensor 
on some fixed super- horizon scale (physical scale L^a(i)). The smoothing scale must be 
below the shortest scale of interest, and to encompass the biggest possible range of scales it 
should not be much bigger than the Hubble scale at the end of infiation. Then we assume, 
by virtue of the smoothing, that spatial gradients have a negligible effect on the evolution 
of the stress-energy tensor. This is the separate universe assumption, that is useful also 
for the evolution of other quantities (we will invoke it later for scalar fields). The separate 
universe assumption will be valid if relevant spatial gradients are negligible, which will be 
the case if is sufficiently far outside the horizon. One generally assumes that the separate 
universes are isotropic (which in our case means that the energy-momentum tensor and the 
expansion are isotropic) so that the separate universes are Robertson- Walker as opposed to 
more general Bianchi universes. That will be the case if scalar field inflation sets the initial 
condition for the Universe, as we are assuming in this paper. The following treatment of C 
goes through however even if the separate universes are anisotropic as might be the case if 
vector fields are involved Q. 

By virtue of the separate universe assumption, the energy continuity equation ^ applies 
at each location. Since we have chosen the slicing of uniform p, the local continuity equation 
reads 

m = -3^|^[p(0+p(x,t)] (19) 

= -3 [H{t) + C(x, i)] [p{t) + p{t) + (5p„ad(x, t)] , (20) 

where ^Pnad is the pressure perturbation on the slicing of uniform density (non-adiabatic 
pressure perturbation). Subtracting the original continuity equation ^ gives 

C(x,t) = ~H{t) ,^ '^Pnad(x^) ^ 

During any era when p{p) is a unique function (the same at each location), ^Pnad vanishes 
and C is time-independent. 

To first order in (5pnad we have, 

C(x,i) = - ,^f_[^\^/ Pnad(x,^), (22) 
p{t)+p{t) 



with 



(5pnad(x, t) = 5p{x, t) - ^'5p(x, t), (23) 

which is valid in any gauge. Also, using Eqs. ^ and ([T^. 

ax,t)^H5t^-HSp/p^\^, (24) 

where 5p is the perturbation on the slicing of uniform a(x, t) (flat slicing). 

A second-order calculation of C is needed only to treat very small non-gaussianity corre- 
sponding to reduced bispectrum |/nl| ^ 1- On cosmological scales, such non-gaussianity 
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will eventually be measurable (and is expected if C comes from a curvaton-type mechanism 
[i,[ni,[il)- But there is no hope of detecting such non-gaussianity on much smaller scales. 

Before continuing we emphasize the following point. A realistic Robertson- Walker universe 
will have inhomogeneities, and one has to take their spatial average in order to determine 
the evolution of the spatially-homogeneous scale factor. This is true in particular for each of 
the separate universes. Various kinds of small-scale inhomogeneity could contribute to the 
spatially-averaged energy density and pressure. The generation of the classical perturbation 
of light fields from the vacuum continues after the smoothing scale leaves the horizon until 
at least the end of inflation, on successively smaller scales. The effect of this perturbation 
on the energy density and pressure during inflation is removed by the smoothing (because it 
is linear to high accuracy in the field perturbation) but the light field perturbations existing 
at the end of inflation might have an effect on the subsequent evolution, for instance on 
massless preheating [l^.'^^ There could be cosmic strings or other localized energy density 
configurations, whose spacing must be much less than the smoothing scale for the separate 
universe assumption to be valid. Finally, through preheating, both heavy and light fields 
can acquire classical perturbations from the vacuum fluctuation on sub-horizon scales. It is 
this effect that will be important for us, as it was for the ordinary preheating calculations 
01 ■ calculate , one must first calculate p and p including the sub- horizon scale field 
perturbations, and then smooth p and p on a super-horizon scale. 



III. HYBRID INFLATION 



In this paper we consider scalar fields with canonical kinetic terms. Their contributions 
to the energy density and pressure are 

P = ^+^E'^' + ^Ei^<^»i' (25) 
p = -^ + ^E'^' + JEi^'^«i'' (26) 

where V is the scalar field potential. Ignoring the metric perturbation (back-reaction) their 
field equations are 

dV 

4 + 3H4 + - = 0. (27) 



A. Slow-roll inflation 



We consider single-field slow-roll inflation, in which the inflaton field 4> has negligible 
interaction with any other field, and is the only one with time-dependence. Consider first 
the unperturbed (spatially homogeneous) universe. The energy density and pressure are 

P0 - V{4') + \4>^ (28) 

= -V{<p) + \4>\, (29) 

and the field equation is 

+ 3i?0 + V'{(t)) = 0, . (30) 
Slow-roll inflation corresponds to almost exponential expansion: 

< ff2 ^ < (slow roll), (31) 



In that case one may have to worry about which reahzation of the ensemble of fluctuations is the observable 
universe [l3[. 
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which leads to the approximations 

ZMIh'^ ^V{(j}), 3H(j) = -V'{(j)) . (32) 

(An over-dot means differentiation with respect to t while a prime means differentiation 
with respect to the displayed argument.) The first derivative of Eq. ([31]) is also supposed to 
be valid, leading to 

101 < H\(j>\, \V"\ < H\ (33) 

For some purposes one or more higher derivatives of Eq. (|3T1) are also supposed to be valid 
leading to further relations, but we shall not need them. 

Ignoring back-reaction, the first-order perturbation Scf) satisfies 

50k + 3-ff(5</.k + {k/af (50k = -V""(0(t))<50k. (34) 

Back-reaction vanishes in the slow-roll limit 0(i) — )■ 0, in any gauge whose slicing is non- 
singular in that limit The slicing of the widely used longitudinal gauge is non-singular, 
and so is the fiat slicing Q. Assuming Einstein gravity, back-reaction is a small effect in 
these gauges. The slicing with uniform energy density is singular in the limit, but Scj) then 
vanishes. 

As each scale leaves the horizon during slow-roll inflation, the vacuum fluctuation gener- 
ates a nearly gaussian classical perturbation (50k with spectrum {H/2tt)'^ . We denote the 
contribution of this perturbation to the curvature perturbation by C,^ (to flrst order in the 
held perturbations, it is the only contribution during single- held inflation). After smoothing 
on a super-horizon scale, Eq. pO[) is valid at each location, with the overdot meaning 
d/dt-pT- and H = [(ia(x, tpr)/(i<:pr] /a (x, tp^), where tp^ is the proper time.*^ 

By virtue of the slow-roll approximation is determined by 0, which means C,^ is con- 
served. It also means that the slicing of uniform p is the same as the slicing of uniform 0, 
so that to flrst order in 50 

C0 = -H5(j,/^. (35) 

This gives 

V^^{k) ^ (h^ I2^4)\ (36) 
where the right hand side is evaluated at horizon exit. 



B. Hybrid inflation 

Hybrid inflation was proposed and named in [l5l | as a way of getting the low inflation 
energy scale that might be required by the axion isocurvature perturbation fiol [l6l | T 7|. 
(See 18| for earlier realizations of hybrid inflation.) It was subsequently found 



be a powerful tool for model building especially in the context of supersymmetry."*^^ Until 
hybrid inflation nears its end a 'waterfall fleld' Xi which has nonzero vev, is flxed at the origin 
by its interaction with 0, up to a vacuum fluctuation which is ignored. This displacement 
of X from its vev gives a constant contribution to V which dominates the total, leading to 
single- fleld slow-roll inflation. After falls below some critical value 0c, the waterfall fleld 
develops a nonzero value x(x, t) which eventually ends inflation. That process is called the 



To first order in 5<f>, the effect of back-reaction on Eq. II34I I takes care of tlie difference between H(t) and 
H{x,t) and between tp^ and coordinate time (see for instance Eq. (5) of [13 ). 

More recently, there is interest in hybrid inflation with a non-canonical kinetic term, in particular DBI 
inflation [2^. The only essential feature, needed to make sense of the hybrid inflation paradigm, is that 
the potential dominates the energy density. 



7 



waterfall; in other words, the waterfall begins when (p falls through its critical value, and 
ends when inflation ends. 
We will adopt the potential 

Vi^, x) ^ Vo + V{<j>) + ^m\<P)x^ + jXx^ (37) 
m\cf>) = g^4>^ -m' = g^ (0^ _ ^l) , (38) 

with < A <C 1 and < g <C 1. The potential is invoked for < (t)ohs^ the value when 
the observable universe leaves the horizon. The effective mass-squared m?{4>) goes negative 
when falls below 0c = w/g, the waterfall then commencing. 

The vev of 4> vanishes and we take V^(0) to vanish at the vev. The inflationary potential 
is supposed to be dominated by Vq which means V{(f)) <C Vq. The requirements that V and 
dV/dx vanish in the vacuum give the vev xo and the inflation scale Vb ~ 3MpH^: 

= Vo = ^^3M^H^ (39) 

The following relations are useful: 

To justify the omission of higher powers of x '^^ '^iH assume xo ^ Mp. This implies 
m/H ^ 1, which is in any case the standard assumption for hybrid inflation. Also, to 
justify the omission of high powers of in 1^(0), we will assume 0obs <C Mp. That requires 

(j)c = m/g^Mp (small 0). (41) 

It also requires [23 | a tensor perturbation well below the present observational limit, implying 

H/Mp < 10'^ (tensor bound). (42) 



Successful BBN requires an inflation scale y/MpH > MeV corresponding to 

H/Mp > IQ-^^ (BBN), (43) 

but viable models of the early universe generally require a far higher value. 

The potential ([571) was proposed in [1^, with V{(j}) = m^0^/2. If one demands ( ~ (^^, 
it gives spectral index n > 1 in contradiction with observation. Many forms of V{4>) have 
been proposed that are consistent with C — S j ^^.d we will not assume any particular 
form. 

Minor variants of Eq. ([37t would make little difference to our analysis. The interaction 
(^202^2 j^jg}j^ replaced by (fp' x^^"^ / where A is a uv cutoff, or the term Xx^ might be 
replaced by ^^^"/A". For our purpose, these variants are equivalent to allowing (respec- 
tively) g and A to be many orders of magnitude below unity. Also, cf) might have two or more 
components that vary during inflation. Most of our treatment of the waterfall will apply 
to that case, if <f> is the field pointing along the inflationary trajectory when the waterfall 
commences. 

More drastic modifications are also possible, including inverted hybrid inflation (26j where 
is increasing during inflation, and mutated/smooth hybrid inflation [8] where the waterfall 
field varies during inflation. Also, the waterfall potential might have a local minimum at 



Taking account of the inhomogeneity in the waterfall begins at each location when </'(x, t) = (pc = raj g. 
Since m ^ H. But the perturbation 5(j> ^ H is small compared with (f>c, which means that this is a small 
effect. 

At least for the usually-considered potentials, this ensures that the waterfall takes at most a few Hubble 
times. In the opposite case the the waterfall can end with an era of two-field slow-roll inflation, involving 
X as well as <p. 
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the origin so that the waterfall proceeds by bubble formation [19|, |27| . Our analysis does 
not apply to those cases. 

In this account of hybrid inflation we have taken and x to be real fields. More generally 
they may correspond to directions in a field space that provides a representation of some 
non-Abelian symmetry group (the GUT symmetry, for the waterfall field of GUT inflation). 
That introduces some numerical factors without changing the structure of the equations. 
For the waterfall field it also avoids the formation of domain walls at locations where x = 0, 
which would be fatal to the cosmology. There may instead be cosmic strings, which are 
harmless if the inflation scale is not too high. For clarity we pretend that 4> and x a-re real 
fields. 



IV. WATERFALL FIELD x 
A. Linear era 

Most papers on the waterfall make two basic assumptions. First, back-reaction is ignored 
so that 

4> + 3H<p-V'^(t) = ~V' {(!))- g'^x'^'l) (44) 
X + Si/x-V^X = -m2(0)x-Ax^ (45) 

Second, the waterfall is assumed to begin with an era when the last term of each equation 
is negligible; during that era, the evolution of x is linear and we call it the the linear era. 
Since </) is smooth on the horizon scale, its negligible evolution at each location is that of 
an unperturbed universe (ie. its spatial gradient is negligible). Choosing a gauge whose 
slicing corresponds to uniform 0, 

Xk + 3ilxk + [{k/af + m\^{t))] Xk = 0. (46) 

Even if slow-roll fails during the waterfall, the slow-roll initial condition ensures that </> 
continues to be determined by 0, which means that the slicing is also one of uniform and 

The energy density and pressure of x are 

Px = "^'Wx' + ^x' + ^IVxP (47) 
= -m\<l,)x' + lx' + l\'^x\^- (48) 

B. Regime of parameter space for our calculation 

In this paper we focus on the simplest regime of parameter space, which will be identified 
in Section [V] In that regime the following conditions are satisfied; (i) slow-roll inflation 
continues until the end of the linear era, (ii) that era takes much less than a Hubble time 
(iii) the change in (j) during that era is negligible. The first two conditions imply that the 
change in (p is also negligible. Under these conditions, the assumptions of the previous 
subsection will be justified in Section IVlII 

Since the changes in (p and (j> are negligible, ni^{(j){t)) decreases linearly with time. Setting 
t = at the beginning of the waterfall and using the dimensionless time t = p,t we write 

m2((/)(t)) = -^^t, = -2g^cj)^ = -2.gm0. (49) 



If the waterfall takes several Hubble times, the scales on which (/> is inhomogeneous will start to come 
inside the horizon but that effect will not be very significant. 
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Our assumptions mean that we can ignore the changes in H and a during the hnear era, 
and we set a = 1. The assumption that the hnear era takes much less than a Hubble time is 

{H|^i)T,,, « 1 [HU « 1). (50) 

The assumption (/> ~ (/)c corresponds to 

(pi/mf « T„-i (0„i ~ 0c). (51) 

We are implicitly assuming that Xk, generated from the vacuum fluctuation, can be treated 
as a classical field, and we shall see that the condition for this is Tni ^ 1. As a result, Eqs. ([50| 
and ([51]) require the hierarchy H <^ ^ <^ m. 

To see when the linear era ends, let us include the small second terms on the right hand 
sides of Eqs. (|44l) and (|45|) . Since the change in is negligible and the perturbation 5(j) is 
small, the unperturbed <j){t) (the spatial average) satisfies 

+ 3iJ<^~3i/</) = -F'(0e)-g'0c(x'), (52) 
where (x^) is the spatial average of x^- Subtracting this from the full equation for (j) we find 

<5> + 3iJ(5'0 - V'<5</) = -l/"((/)c)50 - g^x^5(t> - g^cSx", (53) 

where 

<5x'^x'-(x')- (54) 

Finally, Eq. (|45l) becomes 

X + 3Hx- V^x = -m^(m)x - Ax' - S{m^)x (55) 

where S{m'^) = m'^{ip{x,t)) — m?{(t){t)). The linear era ends when the first term on the 
right hand side ceases to dominate, for one of these equations. For Eq. ([52)) this happens at 
roughly the epoch when (x^) has the value 

Xly = m\^\/gm. (56) 

The constraints (|3ip and ()4ip imply Xni ^ XO; which means that the second term on the 
right hand side of Eq. ([55]) is smaller than the first one at this epoch. We will show in 
Section [VlII that the last terms of Eqs. ([55]) and ([55]) can also be ignored, which means that 
Xni marks the end of the linear era. Afterward, the last term of Eq. ([52[) will cause (/)(x, t) 
to quickly decrease, and m'^{(j)) will quickly move to its final value ending inflation and (by 
definition) the waterfall. 



C. Solution of the waterfall field equation 

By virtue of our condition (ii) we can ignore the expansion of the universe, setting H = 
in Eq. (UHl) to get 

= x{r, fc)xk(r), x^r- = (m^ (r) + fc^) / (57) 

The solutions are the Airy functions Ai(x) and Bi(a;). 

In the quantum theory, Xk becomes an operator Xk. Working in the Heisenberg picture 
we write 

Xk(T) = Xfe('r)ak + Xfc(T)a-k (58) 
[ak,ap] = (27r)353(j^_p) (59) 

The mode function Xk satisfies Eq. ([57[) and its Wronskian is normalized to —1. We choose 

Xk = W2^ [Bi(x) + iAi{x)] . (60) 
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Then, in the early-time regime x —1, = m'^{T) + fc^ is slowly varying {\Ek \ ^ El) and 
the theory describes particles with mass-squared m^it). Indeed, the solution in this regime 
is 

Xfe(r) = (2^)-l/2|x|-V4e--/Vi^l^l^'^ (61) 

which can be written 

Xfe(^)--i=e-/*^*^^ (62) 

We choose the state vector as the vacuum, such that ak| 0| >. A significant occupation 
number is excluded since the resulting positive pressure would spoil inflation [28| . 



D. Classical waterfall field x(x,t) 

We will assume that t becomes much bigger than 1 during the linear era. Then there 
exists a regime x = t — (A://i)^ ^ 1 in which 

Xk{r)c^i2f,)-'/'x-'/^e^^''\ X. ^ Mv^Xfe, (63) 

the errors vanishing in the limit x — > oo. In this regime Xk{T) — Xfc(T)(ak + o-k) to 
high accuracy. As a result, Xk is a practically constant operator times a c-number, which 
means that Xk is a classical quantity in the WKB sense. By this, we mean that a suitable 
measurement of Xk at a given time will give a state that corresponds to a practically definite 
value Xk far into future (We have nothing to say about the cosmic Schrodinger's Cat 
problem that now presents itself.) After the measurement, Xk(''") is classical and 

Xk(T) oc Xfc(T). (64) 

Keeping only the classical modes Xk, we generate a classical field x(xiO- Its spatial 
average vanishes and we can treat it as a gaussian perturbation. Its spectrum is defined 
by Eq. (|5]), with the ensemble corresponding to the different outcomes of the measurement 
that is supposed to have been made to produce Xk(i)- The expectation value in Eq. ([8]) can 
therefore be identified with the vacuum expectation value (xkXp)j which gives P^(fc, r) = 
X^(t). The classical mode function Xk{T) is nonzero only in the regime x = t — / ^ 1, 
where it is given by Eq. (1631) . At fc we have 

PxiO.r) = xLoW = (2Mr^/^)-^e^^'\ (65) 
At fixed r, is a decreasing function of k. At fc^ <C /i^r, 

PAk.r) = X^(r)^xLoWe-'^''/"' (66) 
= xLo(r)e-''^''^^\ fc^(r)^MV2V^, (67) 

which means that Px{k) decreases exponentially in the regime k'i{T) <C fc^ ^ /x^r. Modes 
with k ^ k^, (r) can therefore be ignored. The dominant modes (corresponding to the peak 
of at fixed r) have k ~ fc*(r). From Eq. (|50p we have fc*(T) ^ H which means that the 
dominant modes are sub-horizon. 
From Eq. p2p we see that#i° 

{x'ir))-PAO,r)klir). (68) 



#10 'j'jjg argument t in (x^(''")) is inserted to remind us that the expectation value of (defined as the 
ensemble average or equivalently the spatial average for a given ensemble) depends only on time and not 
on position. 
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Doing the integral we find 

(x'(r)) = ^Px(0,r) dfcfc2e-('=Vfc^M)^(2^)-3/2p^(0^^)fc3(^)_ (69) 

From Eqs. and Tni ^ (ln(xni/M))^^^> with Xni given by Eq. ([55]) . To have a 

classical era r 3> 1 we require Tni ^ 1. To get an upper bound on t^i we use Xni <C Xo ^ Mp 
with n ^ H, to find Tni < (ln(Mp/iJ))^/^. According to Eqs. (gD) and (031), this upper 
bound on Tni is of order 5 to 20 with the lower end of the range far more likely. 

Since the dominant modes have ~ k'^i''') ^ /^^'''i we can set a; ~ r to get Xfc — Mv^Xfe- 
At a typical location we have to a good approximation 

x(x,t) =/xV7x(x,t). (70) 

We also have 

(|VxP) ^ 1 d^kk'\Xk\^ ^ kUr){x') « (x') = f^Mx')- (71) 

At a typical location, | Vx| is of order /c*x, which is negligible compared with the typical 
time-derivative jjLy/rx- Integrating Eq. ([70)) over an extended time interval gives 

x(x,T)cxe(2/3)-'''. (72) 

According to Eqs. (jTO)) . and ([7T|) . this approximation ignores the time-dependence of 

/c* (and less seriously, the prefactor in Eq. ((65|) '). That approximation is good over a time 
interval Ar ^ t, and we will use it over such an interval ending at r ~ Tni. 

These results may not hold near places where x(x, r) vanishes, because there is then a 
cancellation between different Fourier components and the relation Xk — /^V^Xfc cannot be 
expected to give Eq. ((70)) . The places where x vanishes are those in which cosmologically 
disastrous domain walls would form if x were really a single field and where cosmic strings 
might form in a more realistic case. At a typical location Vx/x is of order fc*, therefore the 
typical spacing between these places is of order k~^ . But at a distance L from one of the 

1 /2 

places, X is typically of order iv| Vx| ~ Lk^,{x'^) , which is well below its typical value only 
within a small region L 't^ k~^. Outside this region, there is no strong cancellation between 
Fourier components, and we may expect the approximations to be valid. 



V. CONSTRAINTS ON THE PARAMETER SPACE 



In this section we identify the part of parameter space in which our assumptions are con- 
sistent. There are four independent parameters, which we will take to be the dimensionless 
quantities*^^ 

5 < 1, Hp = H/Mp < 1, H^ = H/m < 1, (73) 

and 

/= (5 X 10-5)"'i/V2^'^= (5 X lO-^j^'T'^'f (fccnd). (74) 

Inflation models are usually constructed so that Vc^^ accounts for the observed Vc^ on 
cosmological scales. Then, if Vc^^ is nearly scale- independent we will have f ^ 1. More 
generally there is an upper bound 

/ < 2 X 10=^ (75) 



Strictly speaking, we cannot consider C0(''^ond)i because C, is usefully defined only after smoothing on 
a super-horizon scale. But the analogue of C, defined on the slicing orthogonal to comoving worldlines 
(usually denoted by TZ) remains useful with a smaller smoothing scale and Eq. 1)7411 as well as the black 
hole bound are valid with replaced by Vf^ [29(1 . 
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corresponding to the bound Vc^^ 10 ^ that is required to avoid excessive black hole 

production from C<P-'^^'^ 

In terms of these four parameters, useful relations are 



Tnl 



ln(i7,„.g-V-i/5)l'^', (76) 



,3 



^ W'gH'H-,'r\ (77) 

and the various constraints become 

10-"^ < Hp< 10^^ (BBN & tensor) (78) 
Hp < IQ-^f (slow roll) (79) 
gHp < H,n (small (j>) (80) 
Hp « Hi (A « 1) (81) 
{Tni/WfH„J « g (i?t„i « 1) (82) 
W\^('gHl « / (^„i ~ ^e) (83) 
g/'/' « i/m (r„i » 1) (84) 

(85) 

The constraint Tni ^ 1 corresponds to the existence of a classical regime, which fails if 
gf^^^ ~ Hm- As seen in Section HV CI there is an absolute bound Tni ^ 20 with t^i < 10 
more likely. Eqs. ([82|l and (|84|) require / ^ (10/rni)'^ which is stronger than Eq. (|75|) . If this 
bound on / is saturated, Eqs. (|82|) and (|84| become g-ff^^ ~ 1 and then Eq. ([83|) becomes 
(7 ^ (10Tj^/^)~^. For smaller / the range for gH^ becomes wider and Eq. ([83|) bounds the 
product gH^. 



VI. PRESSURE AND ENERGY DENSITY OF x 

Taking the spatial gradient to be negligible (ie. excluding places where x is close to zero) , 
Eqs. gZl), dm), and dZQl) give 

Px(x,t) ~ -i^Vx'(x,T) + ix'(x,T)~0 (86) 
P^(x,t) ^ iMVx2(x,r) + ix'(x,r) ^ x'(x,r) ^^Vx^ ^ |m2(^)|2x2. (87) 

The cancellation of the terms in the first expression occurs because we have ignored 
the expansion of the universe when calculating x, so that energy conservation prevents 
moving away from its initially zero value. Instead of taking the expansion into account 
directly, we invoke the local energy continuity equation for p^, which will be valid to a good 
approximation because the gradient of x is negligible. This gives p^(x, t) = — 3-ffp^(x, t). 
Using Eq. ([7^ to integrate p gives* 

3_f^ 3Ht 
(x, t ) = - ^—^Px = - ^T^TjPx • (88) 

Since Ht^i <C 1 and t ^ I, this gives \p^\ <^ \p^\ as is required for consistency. 
To a good approximation the spatially averaged pressure is 

(PxM)=A*W(r)) = (x^) (89) 



^^"^ Since </> and x have independent vacuum fluctuations, is the sum of V^^ and T'^^ , and the black hole 
bound appHes to each of them. 

Eq. ([72} will be a good approximation because the integral is dominated by values of t close to 
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and the pressure perturbation is 

Sp^i^,T)=ti^TSx\^,t)=Sx^ (90) 

where 

Sx' (x, r) ^ x' (X, r) - {x' (r)) , x' (x, r) ^ (x, r) - (x^ (r)) . (91) 

Muhiplying Eqs. (f89|) and ()90p by —iH /2ii^ we arrive at the corresponding expressions 
for p^. 

Our goal is to evahiate the contribution of x to the curvature perturbation which 
depends hnearly on (5p^ and (5p^, and hence on 5x^. Since C is by definition smooth on a 
super-horizon scale we are interested only in 5x^ smoothed on such a scale. Let us determine 
its statistical properties, starting with its spectrum. 

Super-horizon scales satisfy a fortiori k <C fc*(T), which means that we can set k = 
within the integral (1161) . This gives roughly 

Ps^2{T,k) ^ Pl{T,Q)kl{T). (92) 

Doing the integral and using Eq. and going to 7^5^^' '^'^'^ 

VsA^,k)^^{x\T)f[k/k^{T)]\ (93) 

Taking the smoothing scale to be the horizon scale, Eq. gives 

J^-r^J^-^k.i-)lH) «1, (94) 

and for a bigger smoothing scale L we should replace H by L^^ , making the bound tighter. 
We see that after smoothing on a super-horizon scale, is almost homogeneous. 

Using the convolution theorem one finds at A: <C fc* (r) the scale-independent expressions 

Bs^2{r)^P^{T,0)kUr), Ts^r) ^ P^{r,0)kl(T), (95) 

etc. After smoothing Sy^ on the horizon scale, 

{{5x^f) ^ P^kl{T)H\ {{Sx')'), - P*kUr)H^ (96) 

etc.. From these follow the skew, kurtosis etc. of the probability distribution of x^, which 
are seen to be small: 

55x2 - (A:*(T)/i7)-3/2 ^ 1, Ks^2 {k,{T)/H)-^ -^1. (97) 

We conclude that after smoothing on a super-horizon scale, Sx^ is almost gaussian. 

In all of this, we have taken the spatial average of x to vanish. That is is true within an 
indefinitely large region because x is supposed to be generated entirely from the vacuum 
fluctuation. But to make contact with observation we should consider a finite box, that is 
not many orders of magnitude bigger than the region in which the observations are made 
[30| . Denote the average within the box by x we have 

X(x) = X + Xcalc(x) (98) 

where the second term is the one we calculated. This gives 

SX^ = -^Xcalc + 2xXcalc. (99) 

For our results to be valid we need ((<5Xcaic))^ to be much bigger than X^(Xcaic)' when ^xLic 
and Xcaic are smoothed on any scale k~^ much smaller than the box size L. The expected 
value of x^ for a random location of the box is (Xcaic) where now Xcaic i^ smoothed on the 
scale L. We therefore expect 



{kLf > 1. (100) 



We conclude that the spatial average x is negligible, unless we live at a a very untypical 
location. 
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VII. JUSTIFYING OUR ASSUMPTIONS 



A. Neglect of back-reaction 



According to the calculation of x in Section ITVl the expansion of the universe is negligible 
and the dominant modes of x have k ^ k^, ^ H. Therefore, the calculation can be regarded 
as a flat spacetime calculation, formulated at each epoch within a locally inertial frame 
corresponding to a box whose size L satisfies H^^ ^ L <C fc^^, that is at rest with respect 
to the cosmic fluid. Since spatial gradients are negligible, each element of the cosmic fluid 
is free-falling. 

Since the calculation within each box takes place in practically flat spacetime, back- 
reaction will have practically no effect on it. It follows that our calculation is valid globally, 
in the gauge with free-falling threads, and with the time coordinate corresponding to proper 
time along each thread (a synchronous gauge). 



B. Neglected terms in Eqs. (fSS]) and (f55]) 



Let us restore the neglected terms in Eqs. ([5^ and (l55t ^^^ We write S4> = d4>i + 6(t>2 and 
X = Xi -I- X2, where (501 and xi a-re the solutions of Eqs. (p4|) and p6)) . that we considered 
earlier. Then Eqs. (|53p and (I55t give 

5'd^2 + iH5'cj)2-VHct>2 = -[V"m))+g-'xl{t)\H2-9''c^c5x\ (101) 
X2 + iHx2-V^X2 = 2g^cbcXiScb2- (102) 

where Sx^ = Xi ~ iXi) before. In writing these equations, we assumed \S(j)2\ ^ 0c and 
1X2 1 ^ Ixili and we will verify the self-consistency of these assumptions. 

Recall that {5x1)' = '^V^^^xk and (^xf)' = ^P^t^xX- Also, \V"\ H'\, H ^ and 
g^Xni g^H^'^/m'^ ^ It follows that the first term on the right hand side of Eq. ()10ip 
can be ignored, and that its solution is 

502 = -f^Sx^ (103) 

Putting this into Eq. (|102p and remembering that xi oc exp[(2/3)r'^/'^] we get 

X2 = 2^0cXi'502/m'- (104) 

T 

Using |(5xP ~ IxiP < Xni one easily verifies that |i502| ^ 0c and |x2| The neglect 

of back-reaction in Eqs. (jlOip and ()102p is justified because the dominant mode of 502 and 
X2 is sub-horizon, being the dominant mode k^ of xi- 

Since X2 is much less than xi, it has a negligible effect on C^. The effect of 502 is also 
negligible, because its contribution 6p2 — —V'6(p2 to Sp is small compared with that of ^x- 



C. Justifying the assumption of a linear era 



The above discussion shows that the assumption of a linear waterfall era is self-consistent, 
and that such an era will end only at the epoch (x^) = Xni when the last term of Eq. ([52| 



We thank Neil Barnaby for pointing out the need to consider these. 
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becomes significant. What about the possibility that, ahhough self-consistent, the assump- 
tion is actually wrong? For the assumption to hold, it is essential that V' dominates the 
evolution during some initial era of the waterfall. In the opposite case where V is com- 
pletely negligible, both and x will be generated entirely by the potential ^^^(^^ — 4'c)x^- 
This case is discussed in [3l| . Noting that the most rapid growth will be in the direction of 
steepest descent of the potential, they find •\/2(0c ~ <!>) — X ^-^d 

V = Vo-^-^x'- (106) 

This potential leads to a non-linear equation for x(x, i), which has been considered in 
[3^ . To handle it, they give (x^) an initial value coming from its vacuum fluctuation. Of 
course that requires a prescription for handling the uv divergence of (x^) that is different 
from ours which is to simply drop modes that are in the quantum regime. (We comment on 
this issue at the end of Section IXlD 

The authors of [s^l find that x grows like exp(fcquantt) with fcquant ^ 9iti- This growth 
rate may be compared with the growth rate x exp[|(T)'^/^] that we found by assuming 
the existence of an initial linear era, where r =_ut with r ^ 1. The latter growth rate will 



be bigger than the former if gm /xt^/^, and 31 1 suggest that this will be the criterion for 
the linear era to exist. That criterion is satisfied with our assumptions, because Eqs. ((77|) . 
rs]) . and ([M)) imply gm ^ We conclude that the assumption of an initial linear era has 



some justification, though there is as yet no conclusive proof. 

VIII. THE WATERFALL CONTRIBUTION TO C 

A. Calculating the waterfall contribution 

The contribution to C generated at a given epoch during the waterfall is 

Cx(x,r) = C(x,r)-C(x,ri), (107) 

where ti ^ 1 is the epoch when x becomes classical. We are interested in C^(x, Tni). 

Following the procedure of @ used for ordinary preheating, let us calculate by inte- 
grating Since is constant, (5pnad receives no contribution from Sp^. Using Eqs. ((28|) . 
(f29ll . ([87]) and ([88]) and remembering that jp^l < bxL we find 

P 

Spnsid = Spx ~ -rSp^ (108) 
SPx- (109) 



02 + {^2) 

Using this result and invoking Eq. (|72p . this gives 



(x2(r)) + <^2 3(x2(r„i)) + 02 3 (x2(ri)) + <^2 ' 

where 

1 5p^{^,T) H dxH^,r) ixHr)} 



From Eqs. §6^ and §9^ we find 



(111) 



(r(T„l))/r =6i/inl. (112) 

We are assuming Htn\ ^ 1, which suggests 6Htn\ ^ 1 hence (x^('''ni)) ^ 0^- Even if GHtni 
is somewhat bigger than 1, (x^) ^ 0^ until r is close to Tni; indeed, from Eq. (|72|) we have 
(x^) ^ 0^ at the epoch r= given by 

T„i - r= ^ ln[(x^i)/02]/2^i/2 ^ ^^^^ (^^3) 
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It follows that the integral (|110p is dominated by the range Tni — r <C t^i, which means that 
Eq. (|72p will be a good approximation. In either case, the integral of Eq. (jllOp is dominated 
by its upper limit so that to a good approximation 

^ 3(x2(r„0)+<^2 ^ ' 

In the regime 6i?ini ^ 1 we have the following interesting result. Here, Eq. (jlOSp 
becomes 

Spnad = - T^^Px> (115) 

Px 

with l^jinadl ^ I^^Pxl- Therefore is (almost) uniform on the slice of uniform p (practically 
coinciding with the slice of uniform p^). That in turn implies that x^(^j ^) almost uniform 
on the slice of uniform p. 

The spectrum of Cx(^j ''"ni) can be written 



2^ {i + mt^if \H 



n.(fc) = ^ . : Z ^ ■ (116) 



This is the main result of our paper. It holds only in the regime of parameter space for 
which i/<ni ^ 1, and only for Tni ^ 1 (so that x is classical). Also, since ^ is defined only 
after smoothing on a super- horizon scale, it holds only for H ^ k. Requiring say (Htni)~^, 
Tni and H/k all bigger than 2 we see that Vq^ is far below the black hole bound Vq ^ 10~^. 

Taken to apply at say (H/k)^ = 10^, the black hole bound implies on a scale k = H 
leaving the horizon N Hubble times before the end of infiation Vq^ ^ e""^^. Considering 
the shortest cosmological scale, the observed value Vt^lk) ^ 10~^ requires N < 7 which is 
unlikely with any reasonable post-inflationary cosmology even with a low infiation scale. 
We conclude that the black hole constraint on the scale k = H almost certainly makes 
negligible on cosmological scales. 



IX. THE 5N APPROACH 



The calculation of (^^ in the previous section uses cosmological perturbation theory. For 
the contributions of light field perturbations to ^, the 6N approach provides a powerful 
alternative making no reference to cosmological perturbation theory. In this section, we 
recall the SN formula and its application to light fields. Then see how it works for the 
waterfall field. 



A. Initial flat slice 



One may extend Eq. PT)) to describe a generic sequence of fixed-t slices, subject only to 
the condition that the initial slice has unit determinant for gij (called a fiat slice) while the 
final slice still has uniform density. Since the local scale factor is uniform on the initial slice, 
C on the final slice is given by 

C(x, t) = 57V(x, t) = 7V(x, t) - N{t), (117) 

where N{x.,t) is the number of e-folds of expansion at a given location, 

/•t(p) 

N(x,t)= dt'd{yL,t')/a{x,t'), (118) 

Jti 



"f^^^ That regime can hardly exist with our restriction Ht^i ■C 1, but it turns out JSS^ that the result still holds 
in a regime Ht^i ^ 1. 
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and N{t) is the expansion in the unperturbed universe with scale factor a{t). The expansion 
is from the initial flat slice at time ti, to the final uniform-p slice at time t, and we suppress 
the argument ti of N because SN is independent of ti. This formula was given to first order 
in ^] and non- linearly in [35].*^^ 

Going to the proper time ipr at each location, Eq. (|118l) becomes 

7V(x,V)^ if (x,ipdi;,. (119) 

Jtip,(x) 

(Without loss of generality one usually takes iipr to be independent of x.) For each of the 
separate universes, tpr is the cosmic time appearing the Robertson- Walker line element, and 
H is the Hubble parameter given by 3AIpH^ ~ p. Since the final slice has homogeneous p, 
the final value of H is also homogeneous. 

B. Initial uniform-density slice 

The contribution to C that is generated between times ti and t is 

C(x, t) - C(x, h) = SN{x, t,h) = iV(x, t,ti)~ N{t, ti) 
i-Kp) 

N{x,t,ti) = / dt'd{x,t')/a{x,t'), 

where both the initial and final slices have uniform p. 
Going to proper time, 

|.tpr(p,x) 

A/'(x,ipr, Vl) = / ^^(X,ipr)dipr, 
Jtp,.(pi,x) 

with both the initial and final values of H determined by the relation 3MpH^ 

C. 5N for the contribution of light fields 

1. Initial flat slice 

Let us now choose the initial epoch to be a very few Hubble times after the smoothing 
scale of C leaves the horizon during inflation, and denote it by . At this stage the light 
fields are smooth on the scale L.^, because their classical perturbation exists only on scales 
much bigger than the horizon (fc ^ aH) while their quantum fluctuation has been dropped. 

One assumes that the values of one or more of the light fields, at the epoch t*, provide at 
each location the initial condition for the evolution of the separate universe. Then, with N 
defined by Eq. (|118l) . we have in terms of these values 

C(x, t)^SN^ Nit, ^*(x), </);(x), • • •) - (Nit, , • • •)>, (123) 

where (f>* without an argument denote the unperturbed values, and the perturbations 64>* = 
(p* (x) — (p* are defined on a flat slice. 

In practice, C is well-approximated by a few terms of a power-series in the fleld perturba- 
tions: 

c(x, t)^Y. w + +^ E (^)^'^* w)'^'^.* (x) + • • • . (124) 



(120) 
(121) 



(122) 



Considering the time-dependence of C, during inflation, which was also the focus in [34| , the first-order 
formula is also given in and its non-linear generalization is given in |37l. l38l . 
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In this formula, the partial derivatives Ni etc. are evaluated at the unperturbed point in 
field space. 

Taking t to be any time after ^ has stopped varying, we arrive at the quantity constrained 
by observation denoted simply by C(x). Since the observed Q is gaussian to high accuracy, 
(I124D should be dominated on cosmological scales by the first term. This term was given in 
[34| . The power series was given in [39|, where it was used to calculate the bispectrum etc. 
that correspond to non-gaussianity in a couple models. 

If there is single-field slow-roll inflation and the only contribution to Q comes from the 
inflaton 0, the perturbation 50* just represents a shift along the inflaton trajectory which 
does not affect the subsequent relationship p{p). Then C, is time- independent, and the first 
term of Eq. (|124p dominates. To first order in 50*, the time shift is 5t = Jt/f* /</>*, which 
with 5N — —HSt gives Eq. (|35l). If there is multi-field slow- roll inflation with two or more 
fields we can put these fields into Eq. ()124p to arrive at ({x, t) during slow- roll inflation. 
For each of the separate universes, the final epoch should be one of fixed y ~ p, which in 
general will mean different final values for the fields. 

Initial untform-p slice 

So far we have taken N to be defined from an initial fiat slice, leading to Eq. (|124p with t^, 
an epoch soon after the smoothing scale Li^ leaves the horizon. This is definitely the formula 
to use if all of the light field contributions become significant during inflation (so that we 
deal with multi- field inflation). There is another scenario though, where a different approach 
works better. This is the curvaton-type scenario, in which the contribution of some light 
field a starts to become significant only at some epoch ti after inflation is over. For this to 
make sense, the smoothing scale chosen for ( has to be outside the horizon at the epoch 
ti . In principle we should keep the effect on p and p of the modes of a that have entered the 
horizon during the era < t < ti, but such modes are likely to have redshifted away (ie. 
\d'^V/da'^\ < is likely during this era) and they are ignored. Also, the coupling of a to 
other fields is supposed to be negligible during this era. The evolution of a at each location 
is then given by the unperturbed field equation, leading to a relation Xiix*) ^nd 

5x(ti, x) = ^5x. + {Sx*f + ■■■■ (125) 

ax* 2 dxi 

The field cr(x, ii) is supposed to set the initial condition at ti so that 

Cx(x, t) = N^{t, t,)5xiti,x) + iiV^x(i, ti) [5x{ti,x)f + ■■■. (126) 

where a subscript denotes a derivative, N-^ = dN/dx etc.. Eqs. (jl25|) and (|126|) give in 
terms of the nearly-gaussian perturbation Sx* which has the spectrum (i?/27r)^. 

Following [3^, this version of the SN formula is generally used for the curvaton-type 
contribution, usually with the assumption that is negligible so that Cx is the observed 
quantity. With that assumption, there is actually no need to assume slow-roll inflation; all 
one needs is the assumption of inflation with a specified Hubble parameter H{t) [llj . 

D. 5N for the waterfall contribution 

Now we show that the SN formula (|122p provides a direct derivation of Eq. (|114p that 
we used to calculate the waterfall field contribution C,^. Our calculation sets H equal to a 
constant. Also, it begins with a slice of uniform 0, and subsequent slices are separated by 



'^^'^ Considering the time-dependence of ( during inflation, which is also the focus in [33l . the first-order formula 
was also given in |3^ and its non-linear generalization was considered (without a specific application) in 

MM- 
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uniform amounts of proper time, so that (j) remains uniform on these shoes and 5p = 6p^. 
Using these results, Eq. (jl22l) indeed gives to first order in (5p^: 

Cx(x,t„i) - H[6t{^,Ui)-St{^,h)] (127) 



P(tnl) 3 p(tnl) +p(inl) 

where 6t is the time interval from the uniform proper time slicing to the uniform p slicing. 
In the second line we dropped 6t{x.,ti) which as seen in Section IVlIII is negligible. 

One might be concerned that setting H equal to a constant is incompatible with 3MpH^ = 
p. But the non-uniformity of H implied by this relation gives a negligible contribution to 
TV: 

SHdt = i H^-P^dt < (129) 



t, 'iJt, P 2 p(i„i) 

which is much less than Eq. (|128p . justify the assumption that H is practically uniform. 

This use of the i5A^ formula looks very different from the one for light fields, but it can 
be made to resemble the latter using the local evolution equation Eq. ([7^ for x- (Recall 
that this equation is needed if we are to use Eq. (|128p within the framework of the present 
paper, because it is invoked for the calculation of p^.) With this equation, the time ini(x) 
required to give the final energy density a uniform value p(ini) is given by 



P(i„0 = P,(^ni(x)) - ^^^;,2(x,n)e^nV^W. (130) 

This gives iV = (_ff//i)Tni(x) as a function of x^(x, ti), and taking its first order perturbation 
we get 



C,(x,p(t)) = .iV^.'^^(^^'^) 



Sx^{^,h). (131) 



Since its input is the same, this will reproduce Eq. (|128p (hence Eq. ()114|) ) as one can check. 
For instance, in the simplest case that /C0(x, Tni) has a negligible perturbation, corresponding 
to (x^(T-ni)) > 0^ it gives 



^ ipN/Hf^ = const - In [x'(x, n)] , (132) 



so that Eq. gives 



. H SxH^,r,) _ H ^x^(x,r„i) 

''^ 2p^, (X2(ri)> 2p^i (x2(r„i)) ' ^ ' 

which indeed coincides with Eq. (I114p when (x^lTni)) ^ 'P'^- 

We emphasize that Eq. (|131l) . which invokes (5x^(x, ti) smoothed on a super- horizon scale, 
is likely to hold only during the linear era. There is no reason to think that the smoothed 
(5x^ will determine the evolution of ( after that era. Rather, in order to calculate that 
evolution, one should calculate p and p keeping all modes of x and and then smooth them 
on a super-horizon scale. 



X. EARLIER CALCULATIONS 

Several papers [3l|, [12, liol - li^ consider the evolution of the waterfall field, without consid- 
ering its effect on the curvature perturbation. Among them, [ill - lisj invoke the assumptions 
laid down in Section [IV Bl We will not comment on these papers, but turn instead to those 
that consider the curvature perturbation. In doing so, we exclude the non-standard scenario 
\m^ \ <t: H considered in [45] and part of (46ll47j. 
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Several papers |48l45ll | treat the waterfall as two-field inflation with nonzero (x)nr For that 
to be valid, we would need to be in the extreme \m^\ <C H scenario (considered in parts 
of [3 W^) where rn^lt) <^ already when the olDservable universe leaves the horizon. 
The calculation in these papers is therefore incorrect, and all of them except [sij reach the 
incorrect conclusion that (^^ depends linearly on x, with nearly scale invariant. (The 

calculation of [si'l, setting without justification (x) = y/Jx^ reaches a basically correct 
result even though the reasoning is incorrect; it makes quadratic in Sx, with Vi^ oc 
and with very very roughly V(^^ ~ 1 at k — H.) 

Two papers |52l [53 | invoke the function N considered at the end of the previous section. 
They consider hybrid inflation with the potential V{(j)) = ^m^^^. The parameters are 
adjusted to give / = 1, and are chosen so that the waterfall takes a small number of Hubble 
times, not much less than a Hubble time as we are assuming. For the present purpose that 
difference can be ignored, because the qualitative evolution of x and is similar in the two 
cases [s^l- With these parameters, (x^(ini)) ^ 0^, which as we saw at the end of Section 
IVHII means that the slicing of uniform p is also one of uniform x- The latter is assumed by 
the authors, without proof. In 52] N is calculated numerically, while in [s^ it is calculated 
analytically. 

Both papers invoke a quantity, denoted by Syl, which is x smoothed on the scale of the 
horizon at the beginning of the waterfall. In [S^], N is regarded as a function of x rather 
than x^, and is inserted into an expression 

{SxLi^,h)f. (134) 

The subscript attached to {N)l indicates that a smoothing procedure has been applied, 
which is necessary because N itself diverges at xl = 0. This divergence occurs because, 
with their choice of parameters, the energy density at the epoch t^i is less than Vq. Such an 
energy density cannot be achieved if x is fixed at 0. 

This expression for (^-^ differs from (|131|) in several ways. First, the expansion is about 
X = whereas it should be about x^ = (x^)- Second (i5xl)^ (which is smoothed and then 
squared) should be replaced by 6x^ (which is squared and then smoothed). Third, {N)l 
should be replaced by the unsmoothed quantity N. As it happens though, none of these 
errors will affect the rough order of magnitude, because with their parameter choice the 
following hold; (i) the dominant mode of x is of order H so that the smoothed quantity Sxl 
will not be very different from x itself, (ii) ignoring the smoothing, the Fourier components of 
(5xi are the same as those of Sx^ (hi) A^(xi) is approximately linear and close to (A^(xi))-L 
away from the origin. 

In [sl] , the correct Eq. (|13ip is used, except that Sx"^ is replaced by Sxl ■ Because of this 
replacement, P^^2{k) ~ P^(0)fcij is replaced by P^2 ~ P^{0)H^ (since 6xl is smooth on 
the Hubble scale) which means that Vq^ is underestimated by a factor (H/k^,)^, which with 
their parameters is roughly (1/3)'^.^^^ In a footnote of [s^, an argument is given which is 
intended to justify the dropping of short-scale modes when evaluating the perturbation in 
X^. They write X = Xi + XS) and to define xs the universe is divided into boxes with size 
^ H^^. Within each box xs is defined as the usual Fourier Series valid within that box, and 
all Fourier coefficients are generated as usual from the vacuum fluctuation. This proposal is 
not viable, because for a typical realization of the ensemble (presumed to correspond to the 
observed universe) x defined in this way has a large discontinuity at the boundary of each 
box in violation of the field equation. 

In the Appendix of 53| , Cx is calculated by integrating (, but they get a result that is bigger 
than the one from SN, by a factor ln[(x^(tni))/0^(ini)]- The discrepancy arises because they 
identify (5pnad with with Sp^, which makes constant in the regime (x^it)) 3> 0^. As we 
see from Eq. (|109p . Spnad in that regime actually decays exponentially. 



Cx(x) 



dx^ 



We are denoting the dominant scale of their calculated x by A:*, the same symbol that we used for its 
dominant scale in our calculation. 
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The papers we considered so far contain errors of principle. We turn now to papers 
[47I [13]'^^^ where the method is in principle correct, but leads to excessively complicated 
equations. Presumably because of this complication, the expressions for (^-^ in [4^.[5^ are not 
correct. The method is to consider equations that include, in addition to the perturbations 
included in our equations, the second order metric perturbation as well as the the product 
xStj). We demonstrated in Section IVlIl that the latter has a negligible effect and we pointed 
out in Section IIIBI that the former will also have a negligible effect. The motivation for 
including these small extra terms seems to be the following. First, field perturbations are 
defined as perturbations away from zero, so that x is regarded as a perturbation. Second, 
'first order' equations are by definition linear in the field perturbations, so that in this 
order x gives no contribution to p or p. Finally, the next ('second') order is defined to 
involve equations containing terms of the form x^, xy and y^, where a; is a first-order field 
perturbation and y is a first-order metric perturbation. It is the inclusion of the terms 
that causes the main complication in the equations. 

As has been noticed elsewhere [s^ 113, li55], the expression for Cy in [s^ j is manifestly 
incorrect because it is not a local function of x-"^^*^ The expression in |47j (Eq. (57)) is local, 
but still not correct. In particular, it should in the regime 6Htni ^ 1 coincide (to first order 
in Spy,) with Eq. (|114p . whereas it differs by having —m^ instead of our — /x^r = m?{(l)) (and 
by a factor 2). Further calculations using the same expression were done in (STj . 



XI. CONCLUSION 



In this paper we have considered the amount by which C, changes during the linear 
evolution of the waterfall field. We have calculated its spectrum in the regime of parameter 
space in which slow-roll inflation continues during the linear era, with the waterfall mass- 
squared decreasing linearly and the waterfall takes much less Ht^x ^ 1 Hubble times. The 
spectrum, given by Eq. (|116p , is proportional to and amply satisfies the black hole bound 
T-'c, ^ 10~^ on those super-horizon scales k ^ H ior which ( is defined. The full parameter 
space for standard hybrid inflation will be explored elsewhere [s^, using instead of ( the 
quantity TZ mentioned before Eq. (I75|) . 

We made the fundamental assumption, shared by all previous authors except those of 31| , 
that the waterfall does indeed begin with a the linear era. We have justified that assumption 
using the work of [sij , who invoke in turn the calculation of (32| . But the justification is not 
quite complete, because the latter calculation invokes a nonzero value of (x^) at the start of 
the waterfall, coming from the vacuum fluctuation, whereas we have set that value to zero. 

This brings us to an issue of principle, namely the regularization of uv divergences. We 
adopted the simplest procedure, of dropping at a given epoch all Fourier modes (of the 
waterfall and inflaton fields) that are in the quantum as opposed to the classical regime. 
That procedure will violate at some level the energy continuity equation, which follows from 
the complete field equations including all modes. With the assumptions laid down in Section 
IIVB[ the violation will be negligible because the dominant classical modes (with fc ~ fc*) 
appear as soon as r becomes significantly bigger than 1; as a result, the subsequent creation 
of classical modes, which violates the energy continuity equation, can be ignored. But there 
seems to be no guarantee that the violation will be negligible in all cases where a classical 
field is created from the vacuum. 

More seriously, the procedure of dropping the quantum modes will not make sense at all for 
the region of parameter space in which the waterfall field fails to enter the classical regime. 
(We identified one such a region in Section IVl working with our particular assumptions.) In 
that region, the waterfall becomes an essentially quantum phenomenon requiring a formalism 
quite different from the one that we presented. 

The procedure of dropping the quantum regime is crude, but something certainly has to 



#19 'Yhe paper [4^ builds o n [46| and corrects some errors in the latter. 

Since the equations in [541 are erroneous, we will not consider [s^ which uses them to consider a metric 
perturbation related to 
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be done because the energy density, pressure and mean-square fields diverge if the entire 
quantum regime is kept. Indeed, an ultra-violet cutoff at some scale k/a — A much bigger 
than the particle masses gives a constant energy density and pressure p = Zp = A^/IGtt^ 
which violates the energy continuity equation jSSj].'*^^^ Several alternatives to simply drop- 
piiig th e quantum regime have been proposed that can apply to hybrid inflation (for example 
|4ll443l . [60| in addition to the procedure of [s^l that we mentioned before) and it is not clear 
how to choose between them if they give different results. 
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